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Abstract 

The semiclassical approximation of coherent state path integrals is employed 
to study the dynamics of the Jaynes-Cummings model. Decomposing the 
Hilbert space into subspaces of given excitation quanta above the ground 
state, the semiclassical propagator is shown to describe the exact quantum 
dynamics of the model. We also present a semiclassical approximation that 
does not exploit the special properties of the Jaynes-Cummings Hamiltonian 
and can be extended to more general situations. In this approach the con- 
tribution of the dominant semiclassical paths and the relevant fluctuations 
about them are evaluated. This theory leads to an accurate description of 
spontaneous emission going beyond the usual classical field approximation. 

03.65.Sq, 42.50-p, 32.80-t 
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I. INTRODUCTION 



Since the sixties of the last century the Jaynes-Cummings model U is frequently con- 
sidered as a simple model to describe a two-level atomic system interacting with an electro- 
magnetic field in a cavity; for recent reviews see 011 . Apart from its relevance to quantum 
optics, in particular laser theory, this integrable quantum model also allows to test approxi- 
mative methods by comparing them with the exact result. In particular, the "semiclassical" 
theory has attracted considerable attention where the bosonic field mode is represented by 
classical c-numbers while the two- level atomic system is represented as a quantum spin-| 
Q- || . In this approximation the Heisenberg equations of motion |7]] are replaced by linear 
operator equations for the spin variables and an amplitude equation for the electromagnetic 
field which is driven by the expectation values of the spin operators. Taking the expectation 
value of the Heisenberg equations for the spin variables, the optical Bloch equations emerge 
which describe the dynamics of a classical Bloch vector on the two-sphere f||J. It is well 
known that this "semiclassical" theory provides results that are equivalent to a full quantum 
mechanical treatment if the mean number of bosons is very large and fluctuations in the 



boson number can be neglected [TU 



While this conventional semiclassical approach treats the cavity field just classically, we 
attempt at a semiclassical theory treating both the atomic and electromagnetic subsystems 
on an equal footing. Starting from the full quantum model we focus on the most probable 
paths of the system within the path integral representation and relevant fluctuations about 
them. Within the scope of (spin) coherent state path integrals we obtain a semiclassical 
approximation going beyond the classical field approximation. For instance, the approach 
yields an accurate description of spontaneous emission. 

The paper is organized as follows. In Sec. II we first solve the Jaynes-Cummings model 
exactly with spin coherent state path integrals in a subspace with fixed excitation quanta 
above the ground state. Then, in Sec. Ill, we examine a semiclassical description which does 
not rely on these subspaces and can thus be extended to more complicated Hamiltonians. 
With coherent state path integrals the leading order of the propagator is determined by 
solving the Euler-Lagrange equations for the classical path. In Sec. IV we consider contribu- 
tions from fluctuations about the dominant path and show that they lead to a decay of the 
excited two-level system by spontaneous emission. 



II. THE JAYNES-CUMMINGS MODEL 

The Jaynes-Cummings model is characterized by the Hamiltonian 

H = a t a + (l + A)S z + A(aS+ + a t S_), (1) 

where a is the canonical annihilation operator of a bosonic field mode with frequency u and 
S± = S x ± iS y , S z are operators of a spin-- describing two levels of an atomic system with 
energy difference %uj . There are two dimensionless parameters, the detuning A = (u —u)/u 
and the coupling strength X — g/u. We use units with u — 1 and h = 1. It is well known that 
the Jaynes-Cummings model allows apart from H for another time independent operator Q 

N = a) a + S z , (2) 
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which measures the number of excitation quanta in the system. Hence, the time evolution 
operator is of the form 



U(T) = e 



-iHT _ e -iNT e -iCT 



(3) 



where C = H — N. Representing the spin operators in the eigenbasis of S z formed by the 
eigenvectors | |) and | j), t ne first factor in Eq.(|3]) may be written as 



e~ iNT = e~ m 



taT (e"^| T)(Tl+e + ^|!)a I). 



(4) 



Introducing further the eigenkets of a^a, invariant subspaces are distinguished. In particular 
the kets | f n — 1) = | |)|n — 1) and | { n) = | i)\n) span the subspace with N = (n — \). 
In this subspace the time independent operator C generates SU (2) dynamics. This can be 
seen explicitly by introducing the operators 



^=^(|Tn-l)(|n| + ||n)(Tn-l|) 
J y =~(-|Tn-l>Un| + | jnXTn-ll) 
J z =|(itn-l)(tn-l|-i|n)an|), 



(5) 



describing the angular momentum of a spin-|. In terms of these spin operators we have 

C = 2X^J X + AJ Z} (6) 

and we see that in this subspace C gives indeed rise to pure SU (2) dynamics. Accordingly 
the propagator may be worked out exactly by a semiclassical approach with path integrals 
in the spin coherent state representation 



\&<p) = e -W* e -i*Jv\ | n - 1). 



(7) 



Following the lines of pUf , we write the spin coherent propagator clS ct gularized path 
integral 



<tfVV iCT |tfy> = Jim / dfMexp{iS[$(t),(p{t)]}, 

with the action 

S[#(t),<p(fi)]= [ T dt 
Jo 

Here the operator C is represented as 

C(#,<p) = {&<p\C\#<p). 

= \^/nsm( y '^) cos(ip) H cos($). 
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~coB(#)<p-C{0,<p) 



(8) 



(9) 



(10) 



The spherical Wiener measure [12,1 
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exp 



t=o 



4u Jo 



dt 



d 2 + sin 2 (0) ^ 



(11) 



enforces continuous Brownian motion paths on the sphere (M is a normalization factor). 
This measure gives rise to a regularization dependent action 

S v [&(t)Mt)\ = jf dt \x v + sm 2 W0 2 ] + ^cos(^)0 - C(tf,</?)}- (12) 

Now, in the semiclassical expansion, we separate the paths cos($) = cos(i9 c ;) + x/y/s and 
ip = (fid + yj\fs in their classical parts and fluctuations around them. The formal limit 
of large spin s — > oo expresses the classical limit. To lowest order, in the Dominant Path 
Approximation (DOPA), the semiclassical expansion gives 

e iS ^ = expj -i£dtC(ti"(t),<p"(t))}(ti"<p"\ti'<p'), (13) 

where ( , d',(p') and ($",{p"), respectively, describe the starting point and endpoint of the 
classical trajectory (d(t),<p(t)), < t < T. For convenience let us introduce the complex 
variables 



( = tan 



Then, the dominant path is determined by 

C = -a^(l-C 2 )+^AC 
T] = i\y/n(l — rf) — iAt], 



(14) 



(15) 



with boundary conditions ((0) = (' and T](T) = rj". Hence, the endpoint of the classical 
trajectory obeys 



ceo 



20 n C'cos(O n T) + % [AC - sin(fi n T) 
2Q n (' cos{Q n T) - i [Xy/nC + A] sin(fi n T) 



V(T) = if' 
with the Rabi frequency 



A 2 

A 2 n + — . 



In terms of the complex variables (0) we get 

C(T) + V , Al-CW 



c(C(rW) 



1 + C(T)77" 2 1 + C(T) V " 



d 



i— log^i + cV'Wcyr) 



A 



Av^(C + r/') + -(1- CVOJ sin(fi n T) 



(16) 



(17) 



(18) 



Now the integral in Eq.(|I~3"D is readily solved and the propagator in the DOPA takes the 
form 



e lS - = a(T) cos ^ j cos ^ j e^"~^ + b{T) cos f ^ j sin e^'+v') 

-b*(T) sin cos e^"^ + a*(T) sin sin e ~^"^'\ (19) 



a(T) = cos(fi n T) - i^- sin(fi n T) 



where 



b(T) = -t^sm(Q n T). (20) 

As discussed elsewhere |TT| for pure SU{2) dynamics the DOPA is exact and Eqs.(|19D,([2"0"D 
give indeed the exact propagator ||. 

In more general situations, such as for the case without rotating wave approximation 



14| , |15|1 , the system cannot be separated into invariant subspaces. Therefore it would be 



interesting to consider a semiclassical expansion that does not rely on the SU(2) generators 

(!)■ 



III. SEMICLASSICAL DYNAMICS WITH COHERENT STATE PATH 

INTEGRALS 



In order to formulate a general semiclassical theory for a coupled spin boson problem we 
make use of product coherent states 



\$(ppq) 



-i(pS z -itfSy i(pQ-qP) I 



T)|0). 



(21) 



These states are generated by momentum and space translations of the normalized vacuum 
state |0) and SU(2) rotations of the S z eigenstate | |). Again, the semiclassical approxima- 
tion is based on the coherent state path integral representation. We express the propagator 

as 



{$"tp"p"q"\U(t)\#tp V Y> = lim / dti a dn b ex.i>{iS\p(t),q(t),d(t),<p(t)]} 



with the action 

S\p(t),q(t),m,v(t)\ 



T 



dt 



-{pq -pq) + - cos($)0 - H (&, if, p, q) 



(22) 



(23) 



For the Jaynes-Cummings model the Hamiltonian takes the form 



H(-d, ip,p,q) = ({}ippq\H\ , dippq) 

!/ 2 2^ 1 + A 

= 2 (P +q) + — 



(24) 



cos(^) + 



A 



2V2 



sm({})e tlp (q + ip) + sin(-(?)e lip (q — ip) 
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Here, the canonical coherent state path integral is regularized by the flat Wiener measure 
0,0 



dfjL a = M a Y[ dp(t)dq(t) exp 



1 f T 



dt 



•2 i -2 

q +p 



t=0 v 2z/ a Jo 

while the spin paths are again regularized by the spherical Wiener measure 

T ( 1 r T r • 

dfi b = M b TT dcos(tf(i))<fy>(t) exp \ / dt $ 2 + sin 2 ($)v? 2 

t =o ' I Av b Jo 1 

These measures give rise to the regularization dependent action 



(25) 



(26) 



Su a ,u b m)^(t),p(t),q(t)} 



T 



dt t 
o I2v a 



•2 i -2 

q +p 



+ 



$ 2 + sin 2 ($)0 2 
+ ^{pq-pq) + ^cos(i?)0 - H(#,(p,p,q)y 



(27) 



In the semi classical expansion we split the paths p = p c i + x a , q = q c \ + y a , cos($) = 
cos(i? c z) + Xb/y/s and <p = ip c i + y^j \/s in their classical parts and fluctuations around them. 
Restricting ourselves to the DOPA, we obtain the propagator 



sin(tf') sin(tf") ( 1 
\l sin(^) sin(tf") exp {" 2 



q"p" - q"p" + q'p' - q'p' 



xexp<i / (it - cos(^)(^ H — (pq — pq) — H($,ip,p,q) 
[Jo 12 2 



(28) 



While for A = this approximation yields the exact propagator, this property is lost for the 



interacting system. Introducing the complex variables fl6|.[T7 



a = —(q + ip) 



( = tan 



rj = tan 



-vp 



(29) 



Eq.(pHl) may be expressed as 



e lScl = exp 



af + \(3"\ 2 - a(T)(3" - a'/3(0) 



2 

c ^y^ji£dt 



;i + C^(o))(i + c(tv) 
\| (i + CV)(i + CV) 



( d/3 _ a/ j) + - # (a, P, C, t/) 



(30) 
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where the Hamiltonian (E4[) reads 



i/ (a, p, C, r]) = a(J -\ — ; ; ; + A- 



2 1 + Cv 



(31) 



Now, the DOPA propagator is determined by the dominant path obeying the classical equa- 
tions of motion 



a = — i 



a + A 
P = i\p + \ 



V 



1 + 07 

c 



T] = —I 

with the boundary conditions 



1 + 07 

C = ^(l + A)C-A(/?-«C 2 

1 + A)rj - \{a - (3rj A ) 



(32) 



a(0) 



V2 



(q' + ip') 



C(0) = tan 



rj(T) = tan y e^' 



(33) 



The system of differential equations (p2D gives rise to a Hamiltonian vector field. Extend- 
ing results in |T7| , IS|| , one sees that this Hamiltonian dynamics is identical to the classical 
mechanics of a spin on the two-sphere coupled to the phase space degrees of freedom of 
a one dimensional harmonic oscillator. Since the covariant divergence of the Hamiltonian 
vector field vanishes, this dynamical system is conservative and no attractor can occur. The 
coupled differential equations (|32|) with conditions (|33| ) express a nonlinear boundary value 
problem. We can find a solution exploiting the invariance of the action (^) under phase 
transformations 



c- 








v - 


-> V e~ iA 






a - 


-> ae~ lA 


(3- 


-> (3e iA . 



The corresponding integral of motion is 



N{a,P,(,v) = aP + 



ii -Cv 



(34) 



(35) 



2l + ( V 

Therefore the Hamiltonian dynamical system becomes integrable by the theorem of Liouville 



Arnold Hl9| . Particularly, by setting u = (1 — (t])/(1 + (r/), we reduce the system to a 
one-dimensional problem of the form 
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with the cubic potential 



The coefficients read 



where 



+ V(u) = 0, (36) 



V(u) = \ 2 (u 6 + a 2 u 2 + aiu + a ). (37) 



a = -2N + 2^ 
AC 

a 2 = 2N+^- 2 , (38) 



C(a,P,C,v) = X aC+ i 3r] + - 1 -^ 1 . (39) 
V ' M ' S ' /; 1 + C?7 2 l + (r] K ' 



Although the potential V(u) is time independent, the boundary values ( |33"D enforce the 
coefficients in Eq.(|3"8~D to depend on the end time T, and the form of V(u) changes with T. 
Next we set v = u + 02/3 and rewrite Eq.fl3"6|) as 



2 

—v 2 = 4v 3 - g 2 v - g 3 . (40) 



This is just the differential equation solved by the Weierstrass elliptic function [20 
p{-^t] Qi\ #3) with the invariants 

92 = -4 fai - ^2 

4/2 \ 

93 = ~3 \9°^ ~ a / 2 ~ 4a °" ( 41 ) 

In the following we suppress these invariants in the list of arguments of the function p. Now, 
the solution of Eq. ( [36]) becomes 

nit) = -^ + p(a 1 + -4*1 , (42) 



v/2 



where 



_! /a 2 1-C>)(0) A ,,,, 

is determined by the inverse Weierstrass function p~ x . Making use of the solution ([£J), the 
equations of motion lead to elliptic integrals which can be solved in terms of the Weierstrass 
elliptic functions p, ( w and a w [EH]. After some algebra one finds for the field coordinates 
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a(t) 



a 



x 



J v w (A 2 + Ax) ~ 

2a'f3(0)<j w (A 2 -A 1 )_ 

exp{-^C»(^2)t}exp{ 



1/2 



X 



A 

+ 2 



2iV + — - + -4=*) iii) 



1/2 



and 



0(0 = 



1 a w (A 2 -A 1 - j=T) 



2a(T)f3»a w (A 2 + Al + ±T) 
A 



1/2 



X 



e W {^=( w (A 2 )(T - t)} exp{-t (l + j) (T - t)} 

<?w{A 2 



± - M 



1/2 



where 



A 2 = p" 1 (f + 2iv) 
The spin variables are found to read 

1 a w (A 3 + A 1 )a w (A 4 + A x ) 



C(t) = C 



1/2 



X 



C'77(0) a w (A 3 - Ax) a w (A 4 - A ± ) 
exp{^=[£ w (A 3 ) + ( W (A 4 )} t} exp{it} 

1 + f - p{Ax + j-t) a w (A 3 -A ± + j.t) a w (A 4 -A 1 + j-tj 
1 - f + p(A 1 + j-t) a w (A 3 + Ax + j-t) a w (A 4 + A Y + ±t) _ 



and 



v(t) = v" 



1 a w (A 3 + Ax + ±T) a w (A 4 + Ax + ±T) 
C(T)V" ° W {A 3 -A,+ JT) a w (A 4 - A 1 + Jt) 

x exp{ A[ Cw (^ 3 ) + C,(A 4 )](T - 0} exp{i(T - t)} 



1/2 



1/2 



X 



1 + «2 

- 1 T 3 



p(Ax + ^) fr w (A 3 - A 1 + j-t) cr w (A 4 - A, + j-t) 

i-f + pjAi + J) + Ax + JO + Ax + jo 



1/2 



where 



4, = P" 1 (^-1) 



A 4 



-if— 
V 3 



The solutions (f44j)- (f49|) give the dominant path in terms of the known initial values a(0) = a' } 
C(0) = (' and final values (3{T) = (3", rj(T) = rj" and as implicit functions of the unknown 
initial values /3(0), T](0) and final values a(T), ({T). Two of these unknowns have to be 
determined numerically. For instance, from Eqs.( 



and fl47|) we obtain two transzendental 
equations for a(T) and ((T) that can be solved by a root search procedure. Then, the two 
other unknowns can be found from the two constants C and N. 

Having determined the semiclassical trajectory, we may insert the result into Eq.(|50"D and 
determine the DOPA-Propagator. Since this propagator obeys a semiclassical Schrodinger 
equation [see Appendix A] , an alternative representation of the propagator reads 



exp 



i i Q dtH(a(t),(3",at),v")\($" V')(P" q"W I')- 



(50) 



With this representation the DOPA propagator is just determined by the the endpoint of 
the classical path. 

Although the dynamical system fl32|) is conservative, it gives rise to stationary states. 
These are the fix points (a, (3,(n,Vn) — (0,0,0,0) and (a, /3, ps, <Js) = (0,0,0,0), where 
p = l/( and cr = 1/77. These points correspond to the states | j 0) and | J, 0) referred 
to as north pole and south pole, henceforth. For a linear stability analysis we just have to 
linearize the spin terms since the equations of motion (^) are already linear in the oscillator 
variables. Expanding about ((n,Vn) we find 



(51) 



and two invariant subspaces in the variables (5(, (5) and (5r], a) appear. The solution satis- 
fying the boundary conditions (B3|) becomes 









/ 1 + A 


-A 


\ 




d 


P 


= i 


A 


1 







~dt 


5r] 










-1-A A 


5r] 




V oc ) 




1 





-A -I ) 





a(t) 

m 
sat) 

8 V (t) 



1 



cosh(fijyT) 

1 

cosh(Q N T) 
1 

cosh(fijvT) 
1 



{a'e-^'cosh [n N (T-t)] 



ir] » e iu> m (T-t) glu 



Hn N t)} 



II e -iu> m (T-t) 



cosh(n N t) - i('e iuJmt sinh [Q N (T - t)}} 



cosh(f27v^) 
with the frequencies 



{C'e^" 1 ' cosh [Q N {T - t)] - i(3"e 
[rf't 



11 -iu m {T~t) 



sinh(fiArt)| 



iu m (T—t) 



0J„ 



COsh(fiArt) 



A 
~2 



1a e 



-iuJmt 



smh{tt N (T -t)]} (52) 



n 



N 



(53) 



Note that for long times the dominant path converges to the corresponding boundary value 
and no oscillations around the north pole take place anymore. 
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In the same way, we linearize the motion around the south pole. Now invariant subspaces 
appear in the variables (5p, a) and (5a, (3) 





(8p\ 


f-l-A 


-A 






(8 P \ 


d 


a 


-A 


-1 







a 


di 


5a 


1 


1 + A 


A 




5a 




V P J 


V o 


A 


l) 




{ P I 



with the solution 



where 



a(t) 


= a'e 


iUmt cos(fi 5 t) - i—e 


- iUmt sin(fi 5 t) 








Pit) 


= P"e 






- t] sin [Q S (T 


~t)} 




Spit) 


1 


~ iulmt cosiest) -ia'e 


~ iuJmt sin(fi s t) 








5a{t) 


i 


-iu m (T-t) cog ]£i s (r_ 


t)] - ip"e~ iuJm{T - 


- t] sin p s (T 




(55) 








1 A 2 

/ A2 + T' 






(56) 



Here, the dominant path does not converge for long times but keeps on oscillating around 
the south pole. North pole and south pole correspond to the local extrema of the cubic 
potential ( j37[) generated by the coupling of the spin-| to a vacuum field. Whenever the field 
becomes filled with bosons, these fix points bifurcate into limit cycles. 

The presence of stationary states leads to strong deviations of the DOPA propagator from 
the exact result for times large compared to oo^ 1 . In fact, for long times the semi classical 
trajectory approaches the saddle point of the cubic potential and stays there for most of 
the time. For the full quantum problem the state | f 0) is not a steady state, rather it will 
decay by spontaneous emission. In the semiclassical approximation spontaneous emission 
arises from fluctuations about the classical path that are neglected in the DOPA. Hence, to 
obtain useful results also for long times, fluctuations about the north pole need to be taken 
into account. 

IV. FLUCTUATIONS 

The semiclassical expansion of the path integral (|8|) leads to second order contributions 
in terms of Gaussian fluctuation path integrals. Denoting by (x a ,y a ) and ixb,yb) deviations 
from the dominant path variables (p, q) and (cos^), ip), the semiclassical approximation 
takes the form 

Lp" p" q"\UiT)\$ tp' p q') sc = e lScl lim / dfi a dfi b exp \i5 2 S[x a (t),y a (t),x b (t),y b (t)]} , 

(57) 
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with the boundary conditions x a (0) = x a (T) = 0, y a (0) = y a (T) = 0, x b (0) = x b (T) = 0, 
Vb(fy — Vb(T) = 0. Since the canonical Wiener measure (|25| ) is of quadratic form, the 
measure of the fluctuation path integral becomes 



1 1 

d/ia = II 7^dx a {t)dy a {t) exp 

t=0 Z7T 



1 



2v a Jo 



dt 



• 2 i -2 

x a + y a 



(58) 



which is of the same form as the original coherent state path measure. On the other hand, the 
spin measure (|^) is not quadratic, and the dominant path ($(t), (p(t)) cannot be separated 
from the fluctuation variables x b and y b . We have 



d Vb = fl 2S A + 1 dx b (t)dy b (t) expJ --!- f 

t=0 47TS { -10 



l r T 



2v b Jo 
cos^)^ 



dt 



sin 2 (t9) 



- cos(^)0 x b y b + 2 TT ; " to + (V , 

sm (17 J v 2 sin (i>J 



+ sin 2 (tf)z/ 2 

2cos 2 ($)^ 2 + sin(2$)^ 



(59) 



and the regularization of the fluctuation path integral becomes in general time dependent. 
However, when the dominant spin path is strictly independent of time, ($(t), f(t)) = ($ , ip ), 
the measure (|59|) simplifies considerably and we get 



2s + 1 



d/ib = ]J dx b (t)dy b (t) exp<{ - 



t=o 



1 



dt 



2u b Jo L s i n 



+ sin 2 {$ )y, 



Then, after a canonical transformation 



(60) 



x b 



x b 



sin($ ) 
y b = sin(i? )|/ 6) 

the measure (p0|) takes for large s the form of the canonical measure (|5 



dfJ, b = J n 7^dx b (t)dy b (t) exp 



i=0 



l r T 



2v b Jo 



dt 



l2 -J2 
X b + Vb 



(61) 



(62) 



Both measures give rise to the regularization dependent second order variational action 



5 S VatVb [x a (t),y a (t),x b (t), y b {t)} 



dt 



o 



, X a + Va) + 7T- ( X b + Vb) 
LV a lv b 



+ ^( x aVa - iaVa) + ^( x by b ~ x byb) - H (x a , y a , x b , y b ,t) 



(63) 



where the Hamiltonian H Q (t) is determined by the second order contributions of the Hamil- 
tonian H expanded around the dominant path 



H (x a ,y a ,x b ,y b ,t) = ai(t)x 2 + a 2 {t)x a y a + a 3 (t)yl 
+&i(t)5 2 + b 2 (t)x b y b + b 3 (t)yl 
+ci{t)x a x b + c 2 (t)x a y b + c 3 (t)y a x b + c 4 (t)y a y b , 



(64) 
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with the coefficients 



1 d 2 H _ d 2 H 1 d 2 H 



sin 2 (^ ) d 2 H 1 d 2 H 1 d 2 H 

l[t > ~ 2s «9cos(tf) 2 ' 2 ^ ~ !d<pdcQs{tiy 3 ^ ~ 2s sin 2 (tf ) <9<^ 2 

c = sin(i9 ) d 2 ff c = 1 <9 2 ff 
dpdcos(-&)' v /isin(^ ) dpdip' 

_ sin(i9 ) a 2 g _ 1 <9 2 ff 

° 3W dqdcoafiy ° 41 J yfisai{# )dqdip' { ' 



For large s, starting and end points in the coherent state fluctuation path integral 
parameterize states \x a (0) y a (0) Xb(0) yb(0)) and \x a (T) y a (T) Xb{T) yb(T)) which correspond 
to product vacuum states. Propagators leading to stationary saddle points (•&„, <p ) may be 
represented now as 

($ <p p"q"\U(T)\$ oVo p'q') sc = e l ^(00|f/ o (T)|00>, (66) 

with the unitary time evolution operator 

U (T) = T t exp{-t£ dtH (t)Y (67) 



determined by the quadratic Hamiltonian 

H a (t) = ai (t) (Ql - ~) + a 2 (t) (P a Q a + Q a P a ) + 03(f) {P 2 a - \) 

+h(t) (Ql -\) + 6 2 (f) (P b Q b + Q b P b ) + h(t) (P 2 - i) 

+d(f) P a P b + c 2 (t) P a Q fe + c 3 (f) Q„P 6 + c 4 (f) Q a Qb, (68) 

describing two driven coupled oscillators. 

As we have seen in the previous section, for the Jaynes-Cummmings model the north 
pole I I 0) becomes a steady state in the DOPA, and it is essential to take fluctuations about 
this state into account. Unfortunately, the description of the spin degrees of freedom with 
spherical coordinates leads to coordinate singularities. Particularly, the azimuthal angle ip 
is undefined at the poles of the two-sphere. To calculate fluctuations about the north pole 
accurately, we change the coordinate system by a rotation. Since rotations are isometrical 
canonical transformations, the spin path measure ( p6|) stays invariant but the kinematical 
term is not preserved. Instead a phase factor appears which only vanishes if starting and 
endpoint of the spin coordinates are identical. 

Within the DOPA, the probability amplitude to remain at the north pole is just a phase 
factor 

(T 0\U(T)\ T 0)dopa = e iS * = exp{-zi^T}, (69) 

and the north pole becomes a steady state. Taking now Gaussian fluctuations into account 
we have 
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(\Q\U(T)\ T0) sc = e^(00|?7 o (T)|00), (70) 
where the the vacuum amplitude is determined by the time independent Hamiltonian 



(P 2 a +QI-1)- (P 6 2 + Ql~l) + \(P b Qa + QbPa). 



(71) 



For convenience we represent the operators Q a , P a and Qb, Pb by corresponding creation and 
annihilation operators a, and b, 



H D = -! + (! + y)(aa f - b ] b) - y(aa f + b ] b) - t\(ab - aW). 



(72) 



Since aa) — b^b commutes with H Q , we rewrite the time evolution operator in the form 

AT 1 A A 

U (T) = exp{-*— }exp{-*(l + -)cJa t) exp{*(l + -)tfbT} U^T), (73) 



with Ui(T) = exp(-iHiT) and 



rA 



(aa f + 6 f 6) +i\(ab - a ] b ] ) 



(74) 



The operators acJ + fetfe, a& and a^tf span the three dimensional su(l, 1) Lie algebra with 
commutators 



a a 



f + b\ ab 



-lab 



a a 



a&, a' 6' 



r/r/ 



f + b ] b. 



(75) 



For this algebra there is a decomposition into one-dimensional SU(1, 1) transformations 
which holds for the whole group, i.e. for all times P2 ]. We start with the ansatz 



C/i(T) = exp{ / u(T)a t 6 t } exp{z/(T)a&} exp{f (T)(aa f + &*&)}, 
which results in the vacuum amplitude 

A, 



(00|f/ o (T)|00) = exp{-*-T + £(T)} 



(76) 



(77) 



Then, requiring that U\(T) obeys the Schrodinger equation d/dTUi(T) = —iHiUi(T), we 
get the relation 



i-(aa j + b ] b) + A(a f 6 f - ab) = fiaW + ve^abe'^^ 
2 



He' 



'■ aW e uab (aa^ 



+ b ] b)( 



i ab —a o'ftt 



(78) 



where we have made use of the Baker- Campbell-Hausdorff formula. Further, the commuta- 
tion relations (|75|) imply 
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e vab {a<J + tfb)e 



+ b ] b)t 



-/x at fit 



(""''•' ah < ' 



= aa f + 6 f 6 + 2z/a& 
= aa f + 6 f 6 - 2 y ua t 6 t 
,V! a& - /i(aa f + b ] b) + /rW. 



(79) 



Now, Eq.flTSP determines the time rate of change of the functions /i, v and £ by the linear 
equations 











-(i 


-1 






-2/i 



-2/i(l — /xz/) 
-2i/ 




(80) 



which are readily solved with the initial conditions /i(0) = 0, v(Q) = and £(0) = 0. In 
particular, we get for the function £(T) in Eq. (ff6|) 



£(T) = *AT + log 
with the Rabi frequency 

n = 



cos(f2T) — i— sin(fiT) 



(81) 



'A 2 + 



A 2 



Hence, the vacuum amplitude (|77|) becomes 

.AT, 
2~J 



(0 0|C/ o (T)|0 0) = exp{i- 



cos(f2T) — i— sin(f2T) 



and the semiclassical propagator with fluctuations 
(t 0\U(T)\ t 0), 



e ^ 



A 



cos(f2T) — i— sin(fiT) 



(82) 



(83) 



(84) 



includes spontaneous emission leading to an instability of the north pole. Eq.fl8I|) gives the 
exact matrix element of the propagator sandwiched between north pole states. 

When the field is initially and finally not in the vacuum state, the semiclassical propagator 
([57]) is no longer characterized by a fix point path. An evaluation of the fluctuations about 
the semiclassical path would then require numerical methods beyond the scope of this article. 
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APPENDIX A: SEMICLASSICAL SCHRODINGER EQUATION 



Here we derive the semiclassical Schrodinger equation for the DOPA propagator given 
in Eq.fl30|). The time rate of change is readily evaluated, and after an integration by parts 
it may be expressed as 



d_ 



1( da(T,T) ,0/3(0, T) (' 

-p — at 



&n(o,T) 

dT 



dT 'I 



{ dT 

da(t,T) 



dT 



1 + C'V(0,T) 1 + ((T,T) V » 



dt 



(3" + a(T,T) 



dp(t,T) 



dt 



9C(t,T) 
dt 



T 



V"-C(T,T) 



dy(t,T) 
dt 



1 + C(T,7V 



-2iH(a(T,T),(3",C(T,T),ri") 



da{t,T) QU ^ UrT1 ,df3(t,T) , ^r)(t,T) - ((^T)^^ 



dT 



-(3(t,T) - a(t,T)- 



dT 



+ 



l + ((T,T)r,» 



J t=o 



dt 



a{t,T) fdf3(t,T) .dH\ p(t,T) fda{t,T) .dH 



dT 



gt,T) 

dT 



dt 

dy(t,T) 
dt 



da 



dT 



dt 



+ i 



dp 



.dH 



— i- 



(l + ((t,T) V (t,T))* d( 

dH 



9C(t,T) 
dt 



+ l- 



dT \ (1 + C(t, T)rj(t, T)) 2 d V 



(Al) 
(A2) 



Using the classical equations of motions, the integral is found to vanish. Then, we rewrite 
the remaining parts in the form 



d_ 

df^ 



-iH(a(T,T),p"X(T,T),rj") 



1 

-a 



g „{ 0a(T,T) | da(t,T) 



dT 



dt 



+ 



da(t,T) 



dT 



0/3(0, T) + 0/3(t,T) 



+a{T,T) 



dT dT 
( dp(t,T) 



dt 



dp(t,T) 



dC(T,T) , d((t,T) 
dT ' dt 



dT 

9C(t,T) 



T + dT 



j -/3(0,T) 
+ C(T,T) (- 



da(t,T) 



dT 

9y(t,T) 
dt 



dy(t,T) 
dT 



+ 



dy(0,T) &n{t,T) 
dT ' dT 



l + ((T,T)rj" 
V(0,T) 9 -^1 



AS a , 



1 + (' V (0,T) 

where most of the terms on the right hand site vanish. Finally we get 



(A3) 



^e iS - = -iH(a(T),p",C(TW)e 



(A4) 
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Note that the matrix element of the Hamiltonian at the endpoint of the dominant path 
(a(T), /3", C(T), rf') generates the time rate of change of the DOPA propagator and not the 
matrix element of the final state \d" tp" p" q") . For a spin-| coupled to a classical field this 
Schrodinger equation generates the exact quantum mechanics [01]]. 
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